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ABSTRACT 

We perform ideal-magnetohydrodynamic axisymmetric simulations of magnetically 
confined mountains on an accreting neutron star, with masses < 0.12Mq. We consider 
two scenarios, in which the mountain sits atop a hard surface or sinks into a soft, fluid 
base. We find that the ellipticity of the star, due to a mountain grown on a hard 
surface, approaches ~ 2 x 1CP 4 for accreted masses > 1.2 x 10 _3 Mq, and that sinking 
reduces the ellipticity by between 25% and 60%. The consequences for gravitational 
radiation from low-mass x-ray binaries are discussed. 
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1 INTRODUCTION 

The magnetic dipole moment fi of a neutron star is ob- 
served to diminish in the long term as the star accretes 
IITaam fc van den Heuvcl 198 a: Ivan den Heuvel fc Bitzarakil 
1 19951 ). although ersj (|l9971 ) argued that fj, may also be 
a function of parameters other than the accreted mass 
M a . The \l-Ma, correlation has been ascribed to a num- 
ber of physica l mechanisms (|Melatos fc Phinnevl 1200 ll : 
ICumm ing 2005]). First, the magnetic field may be dis- 
sipated in the stellar crust by Ohmic decay, accelerated 
by he ating as the accreted plasm a impacts up o n the 
star iKonar fc Bhattacharval 1 19971: lUrpin et all 1 19981 : 



iBrown fc Bildsteq 

magnetic flux tubes may be dragged from 
conducting core by the outward motion of 
vortices, as the star 



19981 : ICumming et al.1 12004 ). Second, 

the super- 
superfluid 



Ruderman et al.1 



spin s down (ISrinivasan et all 
19981 : IKonar fc Bhattacharval 



1990 



1999 



Konenkov fc Gepperdl200ll ). Third, the magnetic field may 



be screened by accretion-induce d curre nt s within the crust 
jBisnovatvi-Kogan fc Romberg! Il974l ; iBlondin fc Freesel 
Il98d : lLovelace et alj 12005^ 7 In particular, the field may be 
buried under a mountain of accreted plasma channelled 
onto the magnetic poles. When M a is large enough, 
the mountain spreads laterally, tr ansporting the polar 
magneti c flux t owards the equato r dHameurv et al.l 19831: 



2001 
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2002: iPavne fc Melatos 2 004; IZhang 
Payne fc Melatos! 120071 : 1 Vigelius i fc Melatosll2008l . 12009b!) 
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IPavne fc Melato"s1 (|2004l ) computed the unique sequence 
of self-consistent, ideal-magnetohydrodynamic (ideal-MHD) 
equilibria that describes the formation of a polar mountain 
by magnetic burial as a function of M a . They found that 
the accreted mountain is confined by the equatorially com- 
pressed magnetic field, which was unaccounted for in pre- 
vious calculations, and that 10 -5 Mq must be accreted to 
lower fi by 10%. Surprisingly, mountains are stable with re- 
spect to axisymmetric ideal-MHD perturbations; they os- 
cillate globally in a superposition of acoustic and Alfven 
modes but rem ain intact due to magne tic line-tying at the 
stellar surface (|Pavne fc Melatodl2007h . The same equilib- 
ria are susceptible to nonaxisymmetric, Parker-like insta- 
bilities (specifically the gravitationally driven, undular sub- 
mode), but the instability preserves a polar mountain when 
it saturates, despite reduci ng the mass ellipticity by ~ 30% 
| Vigelius fc Melatos! 120081 ). Recently, IVigelius fc Melatos! 
(|2009bl ) considered resistive effects. They found that the 
mountain does not relax appreciably for realistic resistivi- 
ties over the lifetime of a low- or high-mass X-ray binary, 
either by global diffusion, resistive g-mode instabilities, or 
reconnection in the equatorial magnetic belt. The Hall drift, 
which exerts a destabilising influence in is olated neutron 
stars (see, e.g.. iRheinhardt fc Geppertll2002l ). is unlikely to 
be impo rtant in accreting neutron stars due to crustal im- 
purities (|Cumming et alj|2004 ICumming||2005l ). 

The investigations outlined in the previous paragraph 
suffer from two limitations. First, the mountain is assumed 
to rest upon a rigid surface. Under this assumption, the ac- 
creting plasma cannot sink into the stellar crust. This is un- 
realistic. During magnetic burial, frozen-in magnetic flux is 
redistributed slowly within the neutron star by the accreted 
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plasma, as it sinks beneath the surface and spreads laterally. 
IChoudhuri fc Konarl (|2002h showed that the time-scale and 
end state of burial are tied to these slow interior motions. 
Second, the accreted plasma is assumed to satisfy an isother- 
mal equation of state. This is an accurate model only for 
neutron stars with low accretion rates M a < 10 _10 MQyr _1 ; 
the thermodynamics of neutron stars accreting near the Ed- 
dington limit (~ 10 _8 MQyr _1 ) is more complicated, with a 
depth-depen dent adiabatic index ([Brown fc Bildstenlll998l : 
iBrownlfeOQOl ). The equatio n of state affects the grow th rate 
of Parker-like instabilities (|Kosiriski fc Hanasj[200rj) . 

In this paper, we seek to overcome the first limita- 
tion. In section [2] we present a new method of compu- 
tationally simulating the growth of a magnetic mountain 
with M a < O.IMq. In section [3j we compare the struc- 
ture of mountains grown on hard and soft surfaces to eval- 
uate the role of sinking. In section [4] the resulting mass 
quadrupole moment is evaluated as a function of M a for 
hard and soft surfaces. A comparison with the results of 
IChoudhuri fc Konarl (I2002T I. and the implications for gravi- 
tational wave emission from rapidly rotating accretors (e.g. 
low-mass X-ray binaries), are discussed in section [5] 



and it fails to accommodate a fluid interior within its fixed- 

bounda ry framework^] 

In iPavne fc Melatosl (|2007h and IVigelius fc Melatosl 
(|2008l ). Grad-Shafranov equilibria are loaded into ZEUS, 
a multi-purpose, time -d ependent, ideal-MHD solver 



2 GROWING A REALISTICALLY SIZED 
MOUNTAIN BY INJECTION 

In order to investigate how a magnetically confined moun- 
tain sinks into the stellar crust, we need a numerical method 
capable of building a stable mountain, with a realistic 
M a , on top of a fluid base. The approach we t ake b uilds 
upon previous work by IPavne fc Melatosl l|2004l . 120071 ) and 
IVigelius fc Melatosl l|200Sl . 1200981 ). Here, as a service to the 
reader, we briefly recapitulate the physical arguments and 
key results from these previous papers, with references to 
the rele vant sections and e quatio ns. 

In IPavne fc Melatosl (|2004h . axisymmetric magnetic 
mountain equilibria are computed by solving an elliptic 
partial differential equation: the Grad-Shadranov equation 
describing hydromagnetic f o rce b alance in axisymmetric 
geometry [Payne fc Melatosl <|2004l ), section 2.1 and equa- 
tion (12)]. The calculation ensures that the mass-magnetic 
flux distribution dM/dxp is treated self-consistently: the 
final dM/dtp is equal to the initial dM/dtf) together 
with the mass-flux distribution of the accreted matter, 
dM^/dtji, which i s cha racterised by the parameter 7/> a 
[Payne fc Melatosl (|2004 ). section 2.2 and equation (13)]. In 
the limit of small M a , the final equilibrium flux solution 
is characterised by the ratio M a /M c , where the character- 
istic mass M c oc M+R^,B+ i s the accreted mass required 
to halve fi Payne & Melatos (2004), section 3.2 and equa- 
tion (30); IPavne fc Melatosl (|2007l ). section 2.2 and equa- 
tion (3)]. The characteristic mass contains the dependence 
of the equilibrium solution on the parameters of the neu- 
tron star, in particular the magnetic field strength B*. The 
Grad-Shafranov equilibria are computed using an iterative 
numerical solver [Payne fc Melatosl (|2004h . section 3.3]; this 
approach only converg es numerically for l ow ac creted masses 
M a < M c « 1O" 4 M [Payne fc Melatosl (|2004T ). section 3.4], 



|Stone fc Normanl ll992alfbl; lHaves et al.ll2006l), and further 
evolved in axisymmetric (Payne & Melatos 2007) and three- 
dimensional geometries I Vigelius fc Melatosl [2008 ) . The 
characteristic mass M c is used to reduce the length scales of 
the simulated neutron star to circumvent numerical difficul- 
ties and render the simu lations compu tationally tractable 
' Payne fc Melatosl (|2007l ). section 3.3; IVigelius fc Melatosl 
2008), section 2.3 and equation (6), and section 4.6]. 
Two approaches are explored to augmenting the mass of 
a Grad-Shafranov mountain, up to M a < 5.6M C : in the 
first approach, additional matter is injected through the 
outer boundary a l ong t he polar flux tube < ip < ipa, 
[Payne fc Melatosl |2007j), section 4.2]; in the second ap- 
proach, the density of the mountain is uniformly increased 
at every point, while the magnetic field is preserved 
[Payne fc Melatosl (|2007h . section 4.4]. A plausible attempt 
to extend this latter approach to include sinking is outlined 
in appendixfX] ultimately this attempt proved unsuccessful, 
and was abandoned. Instead , the method pre sented in 
this section uses ZEUS-MP (|Haves et all 2006) to build 
magnetic mountain equilibria from scra t ch; thi s approach 
was first proposed in IVigelius fc Melatosl (|2009al ). 



2.1 Outline of the method 

The setup of the simulations presented in this paper is de- 
scribed schematically in Figure [1] Three numerical experi- 
ments are performed: growing onto a hard surface, growing 
onto a soft surface by injecting matter from below, and re- 
peating the latter experiment by injecting matter at some 
altitude. 

To simulate accretion, we inject matter from below, 
through the inner boundary of the simulation at r = r m - m . 
One might expect a realistic simulation of accretion to add 
matter from above, through the outer boundary r = r max . 
The two scenarios are, however, equivalent in ideal MHD; 
the magnetic field is frozen into the fluid, which is thus con- 
strained to move along lines of magnetic flux. Provided that 
the simulation reaches equilibrium, it becomes inconsequen- 
tial, with respect to ideal MHD, through which end of a flux 
tube matter is added. This is because matter cannot cross 
flux surfaces in ideal MHD, so the mass column dM(ip) be- 
tween ip and ip + 8i[> adjusts to reach the same hydrostatic 
radial profile in equilibrium, whether it enters slowly from 
below or falls slowly from above. In the presence of gravity, 
which (in the case of a sinking mountain) induces steep den- 
sity gradients in the fluid base, the results to be presented 
in section 1331 confirm that this situation remains true; two 



1 In IPavne fc M elatos ( 2 0041) . the Grad-Shafranov equation is 
solved subject to Dirichlet and Neumann conditions at fixed 
boundaries. Mathematically, one can formulate a well-posed 
boundary-value problem for the Grad-Shafranov equation in the 
presence of a free boundary, e.g. the sinking base of a mountain; in 
practice, however, this is an extremely difficult problem to solve. 



Sinking of a magnetically confined mountain on an accreting neutron star 3 




r„,i n R, r m « r min R, r„ r mi „ = fi', r . 

(a) (b) (c) 

Figure 1. Diagrams illustrating schematically three mountain growth scenarios. The simulation region is bounded by r m ; n < r < r r 



< (9 < 7r/2 (thick lines), and represents a quadrant of the star. Boundary conditions assume symmetry about 8 = and reflection at 
9 = n/2. The surface of the star is located at r = -R*. Three sub-regions are identified. The fluid interior beneath the surface, containing 
mass Afbasci is shaded gray. The region where the mountain mass M a is injected into the simulation is shaded in a gray to white gradient; 
the gray is proportional to the injected flux [see equation Q] as a function of 6. The outer atmosphere of the star, containing mass 
Aiatm, is unshaded. A central gravitational point source is labelled with its mass M*. For mountains grown on a fluid base, M a can be 
injected | (a) | at the inner boundary r = r m ; n , or |(b)| at the stellar surface r = i?*. For mountains grown on a hard surfa ce, |(c)| there is no 
fluid interior; the inner boundary is identical to the stellar surface r = r m i n = JJ* (Afbase = 0). See the text in section \2. 3 1 



different injection scenarios (described below) give elliptici- 
ties consistent to within 10%. There remains, however, the 
subtle and difficult question of irreversible magnetic recon- 
nection at the grid corners, which remains unresolved (see 
the discussion in section [3T4]| . 

In practice, it is advantageous to add matter through 
the inner boundary, because we wish to inject along particu- 
lar flux tubes, and this is easiest to do at r — r m i n , where the 
magnetic footprints are fixed in place (unlike at r = r max ). 
This constraint, known as mag netic line tying, contributes 
to the stability of the mo untain l|Goedbloed fe Poedts|[20o3 ; 
IVigelius fc Melatosl [2008). It is well justified physically, pro- 
vided that r m i n lies deep enough within the star, so that the 
fluid base (and frozen-in magnetic flux) remains relatively 
stationary, and is not significantly perturbed by the spread- 
ing and sinking of the mountain. This is the case if the mass 
Afbase of the fluid base, initially in the region r m i n < r < i?*, 
is much greater than M a . To confirm that the mountain does 
not greatly push the crustal material, we first calculate the 
fraction of Mb asc contained in each grid cell, and then de- 
termine the change in this quantity between the initial and 
final times of the simulation; this gives the change in the 
spatial distribution of Afbase over the simulation, as a func- 
tion of the grid cell. For all simulations with sinking, the 
median change in Mb aso , over all grid cells, is on average 
~ 10%; thus, the distribution of the fluid base does not 
change much during accretion. Recent molecular dynamics 
simulations of crystalline n eutron matter, which pre dict a 
high breaking strain ~ 0.1 ([Horowitz fc Kadaul 12005^ 1. also 
support the line-tying hypothesis. 

When a mountain is grown onto a fluid base Mb asc , 
a difficulty arises. Z EUS-MP models a single fluid, with a 
unique velocity field l)Haves et al.l 1200*61 ); there is no facility 
for simulating the movement of one fluid, the injected moun- 



tain, with respect to another fluid, the stationary crust 
We are left with two alternatives: to assign the same veloc- 
ity to the injected mountain and the crust (the behaviour 
of ZEUS-MP 's "inflow" boundary condition), or to assign a 
negligible or zero velocity to the injected mountain, in order 
to keep the base stationary. In the first case, ZEUS-MP fails 
catastrophically for desirable values of the injection veloc- 
ity (> 5% of the escape velocity). In the second case, which 
we study in section [3] mountains remain subterranean and 
never rise to the stellar surface r — As a check, therefore, 
we examine two scenarios: injection at r = r m i n and r — 7?*. 
We show in section [3T31 that the results in both scenarios are 
quantitatively alike, confirming their robustness. 

Throughout this paper, we adopt the viewpoint that the 
accreted matter and the mountain are one and the same; the 
accreted mass and the mass of the mountain are identical 
and are both denoted by M a . This is a matter of termi- 
nology, not physics. There is no "hard edge" to the moun- 
tain; matter is accreted on all flux surfaces < ip < ip* 
[see equation Q in section [2.3] . not just on the polar cap 
< i(> < ipn, which contains ~ 63% of M a . Under the as- 
sumption of ideal MHD, matter cannot spread across flux 
surfaces, i.e. there is no Ohmic diffusion. We also do not 
model the accreted matter once it has sunk beyond the crust, 
as do e.g IChoudhuri fc Konarl (|2002T ); see the discussion in 
section [5] 



2.2 Initial setup 

The initial setu p of our simulations closely follow s 
iPavne fc Melatosl \200H) and IVigelius fc Melatosl (|2008h . 



2 ZEUS-MP can track the concentrations of comoving compo- 
nents within the same fluid; we exploit this in sectioii [2.3l 
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The setup of ZEUS-MF^ is through a set of parameters 
which control: the geometry of the problem, the physical 
phenomena to be modelled (e.g. MHD, gravity), the sim- 
ulation grid and its boundary conditions, the equation of 
state, and the choice of tim estep. Appropriat e valu es for 
these parameters are given in IPavne fc Melatosl (I2007I). sec- 
tion 3 and appendix Al, and in lVigelius fc Melatosl (|2008). 
sections 2.2-2.3 and appendix A. 

To avoid numerical difficulties with steep magnetic 
field gradients, we simulate a scaled-down neutron star, 
where the mass M* and radius /?* are artificially reduced, 
while the hydrostatic scale heigh t ho = c^R^/GM^ is kept 
constant |Pavne fc M elatos 2 0071 ). The down-scaling trans- 
formation preserves the equilibriu m shape of the moun- 
tain exactly in the small- Af a limit (|Pavne fc Melatosl |2004| . 
l2007h and has been validated a pproximately for M a < 
20M C jVigelius fc Melatosl 120081 ). We use dimensionless 
units within ZEUS-MP, setting the isothermal sound c a and 
gravitational constant G to unity, and adopting ho as the 
unit of length. Table [1] explains how to convert between an 
astrophysical neutron star, the scaled-down model, and di- 
mensionless ZEUS-MP units. 

The simulations are performed on an axisymmetric rect- 
angular grid with N r cells spaced logarithmically in r, and 
Ne — 64 cells spaced linearly in 9. The logarithmic spacing in 
r is determined by the ratio Arjv,.-i/Aro of the maximum to 
minimum radial grid spacing (see appendix[B]). This ratio is 
chosen large enough to concentrate grid resolution near the 
inner boundary, but small enough to ensure reasonable run 
times. We set r max = 1.27?* = 60/io to give the mountain 
ample room to expand without meeting the outer bound- 
ary, and stipulate reflecting boundary conditions at 9 — 
and 9 = 7r/2, "inflow" boundary conditions at r — r m i n , 
and "outflow" boun dary conditions a t r = r max ; more de- 
tails can be found in IPavne fc Melatosl l|2007h . The magnetic 
field is initially that of a dipole, and B* is its magnitude at 
the polar surface. 

A gravitational point source M* is placed at r = 0, and 
self-gravity is ignored. The density field is initialised to be 
the static atmosphere of an isothermal fluid with no self- 
gravity: 



p(t — 0, r) — p Rt exp 



1 

~R~, 



(1) 



Ideally the region r > R+ should start evacuated, but ZEUS- 
MP requires the density to be nonzero everywhere, so we set 
A/atm = 5 x 10 _6 Mq (approximately 4% of the mass of the 
smallest mountain; see Table [2} ■ Integrating equation {TJ 
over the region r > i?* (see Figure [TJ fixes the density at 
the stellar surface pn t in terms of M atm - In contrast, we 
require the mass of the fluid base Mb asc (when a soft surface 
is being modelled) to be much larger than the mass of the 
mountain, as discussed in section [2~T1 In all runs, we choose 
Mbase/Mn. ~ 10. Integrating equation ([T]) over the region 
containing Mb ase then fixes r m i n . 
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Figure 2. (left) The accretion rate M a (t), given by equation J3J. 
(right) The initial mass-flux distribution dM/dip, given by equa- 
tion igj. Contours are at (right to left) 0.05, 0.1, 0.3, 0.5, 0.7, and 
0.9 of the maximum. 



2.3 Injection procedure 

ZEUS-MP's "inflow" boundary condition permits injection 
at the edge of the simulation volume. To enable injection 
at r = R t , as in Figure 1(b) we implement a more flexible 
custom procedure, and use the built-in "inflow" condition 
only to tie the magnetic flux at r = r m i n . We describe the 



procedure briefly below; further details are in appendix [Cj 
We inject mass M a into an injection region ri n j < r < 
rinj+Srinj, < 9 < n/2, over a time interval < t < T a . (We 
set Sri n j = O.l/io throughout.) The flux of accreted matter 
at time t entering a point (r, 9) in the injection region is 
given by 

9 3 M a 



(t,r,9)ocM4t)^(r,9) 



dtdrd9 
where we choose 
M a (t) (xt 2 {T a -t) 2 



and 
dM a 

dip 



(r, 9) oc exp(— bR+r 1 sin 2 9) . 



(2) 



(3) 



(4) 



Version 2.1.2, available 

http : //lca.ucsd. edu/portal/codes/zeusmp2 



from 



The normalisation of equation ((2} is chosen so that, for each 
simulation, the mass of the mountain is equal to M a at time 
t = T a , i.e. M a (t = Ta) = M a . After time t = T a , no further 
mass is added, but we evolve the system up to t = £ max = 
1.5T a to test the stability of the mountain obtained. 

Equation (|3J) determines the rate of accretion; it is plot- 
ted in Figure [2] (left). The functional form was chosen to en- 
sure numerical stability in ZEUS-MP, and has no particular 
astrophysical justification, except to ensure that a mountain 
builds up to its target mass smoothly over the time scale T a . 
For this reason, it is a smooth bell-shaped function, designed 
to avoid any discontinuity in the accretion rate, which might 
excite undesired oscillations in the fluid or provoke numeri- 
cal instabilities. 

Equation (|4| gives a mass - flux d istribution consistent 
with that of IPavne fc Melatosl (|2004l ): it is plotted in Fig- 
ure^ (right). It does not attempt to model the interaction of 
the accreted matter with the magnetosphere, from which the 
mass-flux distribution would originate; instead, it is chosen 
such that the majority (~ 63%) of the accreted matter falls 
on the polar cap < ip < f/)a. The parameter b = i/ 1 */ 1 /^ = 3 
determines the polar cap radius R* sin -1 (6~ 1//2 ). It is deter- 
mined astrophysically by disk-magnetosphere force balance, 
and is related to the ste llar magnetic field via b oc ip+ oc 
l|Pavne fc Melatosll2004l ). In this theoretical paper, however, 
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Table 1. Conversion of physical quantities into dimcnsionlcss variables in the simulations. Physical quantities are first converted to their 
values in the scaled-down model by multiplying by [a(Simulation)/a(Astrophysical neutron star)]™, where a = R*/ho paramctcriscs the 
curvature down-scaling, and n is listed in column 4. Scaled-down physical quantities are then reexpressed in the di mensionless units 
of ZE US-MP according to column 6. The table is divided into three horizontal parts containing: stellar parameters (Payne & M elatosl 
2004), simulation control parameters (see Figure [TJ, and simulation outputs. 



Quantity 


Symbol 


Astrophysical 


Down-scaling „. , ,. 

° simulation 


„„ TT „ .... 
ZEUS-MP dimensionless unit 


neutron star 


index n 






scaling ratio 


a 


1.9 X 10 4 


none 


50 


none 


stellar mass 


M* 


1.4M 


2 


1O _5 M 


M = c^ha/G = 8.1 x 10 24 g 


stellar radius 


H, 


10 6 cm 


1 


2.7 X 10 3 cm 


h = clRl/GM± = 54 cm 


stellar magnetic field 


Bi, 


10 12 G 


none 




B = c^/G 1 / 2 ^ = 7.2 x 10 17 G 


isothermal sound speed 


Cs 


10 s cm s _1 


none 




Cs 


critical mass 


M c 


1.2 x 1O" 4 M 


4 


6.1 x 1O" 15 M 


M 


inner boundary 


^ min 




1 


see Table |2] 


ho 


outer boundary 


''"max 




1 


1.2i?,* 


ho 


accreted mass 






4 


see Table [2] 


Mo 


mass of outer atmosphere 


M atm 




4 


5 x 1O- 6 M 


Mo 


mass of fluid base 


-Afbasc 




4 


10M a 


Mo 


mountain density 


pX a 




1 




po = M /hi = 5.2 X 10 19 g cm" 3 


magnetic field 


B 




none 




Bo 


cllipticity 


e 




2 




none 


time 


t 




none 




to = ho/c s = 5.4 X 10~ 7 s 



we treat b (and therefore B*) as a free parameter, and do not 
attempt a self-co nsistent solution of the disk-magnetosphere 
system (see, e.g.. lRomanova et al.ll2008h . With this freedom, 
b is chosen unrealistically la rge to preserve numerical stabil- 
ity (|Pavne fc Melatosll2004 ). 

We use ZEU S-MP 's multi-species tracking facility 
l|Haves et al.ll2006h to record, throughout the simulation, the 
fraction of the density, < X a (t, r, 6) < 1, that originates 
from accretion (i.e. added at t > via the injection proce- 
dure), as opposed to from the initial configuration at t = 0. 
This allows us to track the spread of the mountain inde- 
pendently of the motion of the remaining (displaced) stellar 
matter. 




Figure 3. Radial component of the magnetic field B r along 6 = 0: 
for simulation S(i?,*,10 2 ) (left), at t = 0.45T a (solid) and t = 
0.85T a (dotted); and for simulation <S(H*,10 3 ) (right), at t = 
0.24T a (solid), and its time of failure t = 0.35T a (dashed). 



We require that the mountain grows quasistatically, in 
the sense that the accretion timescale T a is always much 
greater than tAifven, the characteristic pole-equator cross- 
ing ti me of an Alfven wave. Following IVigelius fc Melatosl 
(120081 ), we compute the crossing time at t — 0, r = 
from the Alfven speed WAifvon = 

(see Table 1), we obtain tAifvcn = irR* /(2«AuV<fa) ~ 400£o- 
The condition tAifven <C T a is verified by comparison with 
the values for T a listed in Table [2] The condition also im- 
plies that the magnetostatic limit always holds: the ratio 
B^M a /A/ a > 7r 3/2 i? + p^ 2 « 2 x 10 8 G s, and from Tables 1 
and 2, B*M a /M a w B*T a > 3 x 10 9 G s. 

For mountains grown on a hard surface, we addition- 
ally set the velocity v(t, r, 6) within the injection region, 
such that the accreted matter is always given a fixed speed 
^inj = 10 _4 c a parallel to the magnetic field B(t, r, 6). The 
value of Winj should be a small fraction of the escape speed 
v csc ~ 4.1c s , so there is negligible mass lost through the 

We find that setting v 



outer boundary (see section 13.2 
carefully is critical to stability. 



3 COMPARING MOUNTAINS GROWN ON 
HARD AND SOFT BASES 

Table [2] lists the parameters of the simulations presented 
in this paper. Mountains grown on a hard surface are la- 
belled 7i(M a /M c ). Mountains grown on a fluid base are la- 
belled <S(ri n j, Ma/Mc), where the injection radius ri n j may 
be either r m j n or R+. The parameters of each run are cho- 
sen to grow a mountain with a particular target mass, 
Mb,. We choose four values for M a in the range 10 -4 - 
10 _ Mq. These values are chosen to demonstrate the ability 
of the injection procedure to generate stable mountains over 
a wide range of masses. This range also encompasses the 
range of M a of real accreting neutron stars (see Table |3}. 
The main source of uncertainty i s the accretion efficiency 
(Ivan den Heuvel fc B itzaraki 1995), which may be as low as 
~ 5% |Tauris et all 2000): this is reflected in the chosen 
range of M a . 

The CPU time required for each run was, on average, 
~ 10~ 2 seconds per grid cell per unit to of simulation time. 
For iS(-R*)-type simulations, one must scale T a with M a to 
prevent numerical instabilities. Even so, run S(R*, 10 3 ) does 
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Table 2. Simulations of magnetic mountains presented in this paper. The accompanying parameters are: the target accreted mass M a , 
in units of M c and Mq; the number of grid cells in the r direction N r ; the radius of the inner radial boundary r m ] n ; the ratio of the 
maximum to minimum radial grid spacing Ar jy r _ i / Aro ; the injection radius r; n j; the injection velocity fj n j; the injection time T a ; and 
the total (successfully completed) simulation time t max . 
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0.35 



not complete; ZEUS-MP aborts at t Ri 0.35T a , when the 
adaptive time-step shrinks below its allowed minimum. Fig- 
ure [3] shows that this behaviour arises when B r diverges at 
r < n n j along the boundary 9 = 0: for Af a = 10 2 M C , B r 
threatens to break out for t < 0.5T a but ultimately set- 
tles down to the equilibrium configuration before t = T a , 
whereas for M a = 10 3 M C , it grows uncontrollably up to the 
time of failure. 



3.1 Verification 

We first check that, for each mountain, (i) we accumulate the 
correct total mass M a , with minimal loss through the outer 
boundary; (ii) the mass above the surface, M atm , remains 
much smaller than M a ; and (iii) for mountains with sinking, 
the mass in the fluid base, Mb ase , remains large compared 
to M a , so that the magnetic line- tying condition at r — r m i„ 



is a good approximation. Figures 4(a) 4(c) and 4(e) show 
Ma(i) = JydVpXa integrated over the simulation volume 
V at time t. We see that the mountains achieve their tar- 
get mass, which remains in the simulation for t > T a . The 
injected mass A/ a (t)/M a , found by integrating equation ([3]) 
with respect to time, is plotted alongside in grey; the two 
curves overlap. Figure 4(b) shows M atm for the hard-surface 
experiment; it is always small. Figures [4(d)1 and [4(f)] show 
Afb a se for the soft-surface experiments; it always exceeds 
x, 10M a , as desired. For all simulations where Mb asc > 0, 
the total fraction of Mb ase l ost through the outer boun dary 
is < 0.01%, consistent with I Vigelius fc Melatos] (|2008l ). 

We next check that accretion takes place in the mag- 
netostatic limit, i.e. that the total kinetic energy = 
J v dVp\v\ 2 /2 is small compared to the total magnetic en- 
ergy E m = J Y dV|B| 2 /(87r). Figure [5] shows the ratio of E k 
to Em as a function of time. We see that E^/Em tends to in- 
crease with M a but typically never rises above 1%, except in 
H(10 3 ) and the incomplete run 5(r m i n , 10 3 ). After accretion 
stops, at t = T a , E^/Em typically falls to less than 10 -4 . 

Magnetic fie ld transport in ZE US-MP is divergence-free 
by construction (|Haves et al.| [2006), but it is worth checking 
whether this property is preserved by the injection algo- 
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Figure 4. Accreted mass M a (t) and masses in the stellar at- 
mosphere M atm (t) and fluid base Mb ase M! plotted in black, for 
M a /M c = 1 (dotted), 10 (dot-dashed), 10 2 (dashed), and 10 3 
(solid). The injected mass from equation J3J is over-plotted in 
gray. The labels beneath each panel indicate a hard-surface (H) 
or soft-surface [S(r m i n ) or 5(i?*)] run; see Table [2] The short 
solid line at the bottom of Figure |(f )| is from the aborted run 
5(R*,10 3 ) 
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(b) 5(r min ) (c) S(R*) 

Figure 5. Kinetic energy E-^(t), normalised by the magnetic 
energy E m (t), for M" a /M c = 1 (dotted), 10 (dot-dashed), 10 2 
(dashed), and 10 3 (solid). The labels beneath each panel indicate 
a hard-surface (7i) or soft-surface (<S(r m j n )) run; see Table [2] 
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Figure 9. Shaded contours of the absolute, normalised difference 
in total density p between (Sfrmin, 10) and S(R *, 10), ove rlaid 
with the pX a (accreted only) contours of Figures [8(b)1 and [8(c)1 
(solid and dashed respectively). The minimum and maximum 
density differences are indicated with small arrows and labelled 
in bold. 




6 

Figure 6. Contours of the absolute, normalised divergence of B 
of mountain ?f(10 3 ) at t = T a . Contour levels are at 5 X 10~ 2 
(solid), 10~ 2 (dashed), 5 X 10~ 3 (dotted). 



placed from its initial configuration, but remains anchored 
to the inner boundary at r = r m i n , demonstrating magnetic 
line tying. We see, in the i/j/"^* ~ 1-2 x 10~ 2 contour, the 
early for mation of the ma g netic "tutu " configuration, ob- 
served in lPavne fc Melatos! (|2004 120071 ) for M a = 1O" 5 M . 

At t = T a , the mountain reaches its target mass (M a = 
10 2 M C in Fi gure 0. Despite sliding towards the equator, 
the accreted matter still exhibits a noticeable variation in 
density with respect to 0; a polar mountain is formed. The 
tutu configuration of the magnetic field is clearly visible; 
see fo r comparison Figures 2 and 4(a) of IPavne fc Melatosl 
(2004). This equilibrium state remains largely unchanged 
when we run the simulation for an additional 0.5T a , during 
which no further mass is added. 



rithm. We find that the mean value of |V ■ B|/ ^QB^/dx') 
is initially < 6 x 10 -3 , and increases by a factor of 3.5 at 
most over the run. Figure[6]shows contours of the normalised 
IV ■ B| for an illustrative mountain. 



3.2 Illustrative example 

We choose mountain 7i(10 2 ), grown on a hard surface with 
M a = 10 2 M C , to illustrate the general evolution of a mag- 
netic mountain during accretion. The top row of Figure [7] 
shows contours of the mountain density pX a , normalised 
by the initial surface density pn t , at four different times. 
The lower row shows the magnetic flux ip, normalised by 
ipi, = B+Rt/2, at the same times. Matter is added predom- 
inately at the pole, as determined by dM/dif). In the early 
stages of accretion (t — 0.1T a ), the magnetic field is only 
slightly disturbed. As accretion progresses, the mountain 
spreads towards the equator, dragging the frozen-in mag- 
netic field with it. The angular span of the tp contours is com- 
pressed from ~ 70° [Figure ffjej] to ~ 20° [Figure fUfj\ . At 
the half-way point (t = 0.5T a ), the flux is significantly dis- 



3.3 Sinking 

We next compare the hard-surface equilibrium state, illus- 
trated in Figures 7(c) and 7(g) with the two experiments 



where we include sinking. Figure [8] shows contours of pX a 
and tp in each of the three scenarios, with M a = 10Af c . The 
hard-surface mountain [Figure 8(a) spreads appreciably, 



and the magnetic flux [Figure 8(d) is significantly displaced 
towards the equator. The density contour pX a //9_R t = 10 
begins a distance ~ 4/io = 215 cm (see Table [TJ above the 
injection radius r — r m i n = 7?* at the pole and sinks below 
the equator to ~ 75% of the polar height of the mountain. 
In contrast, the same contour of the sinking mountain grown 
at r = r m i n [Figure 8(b) begins ~ 0.6ho above the injection 
radius 
just 



r m i n at the pole and sinks below the equator to 
19% of the polar height of the mountain. From the 
pXj/pij, = 10~ 3 contour, we see that the accreted matter 
is confined tor - r m i n < 0.7/io above the inner boundary at 
the pole and r — r m i n < 0.15/io at the equator. 

In short, the sunk mountain grown at r = r m i n hugs 
the inner boundary and pole and resembles the initial mass 
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Figure 7. Hydromagnetic structure of hard-surface mountain 7i(10 2 ) at times t/T a = 0.1, 0.5, 1.0, and 1.5. (Top row) Contours of 
accreted density pX a / pn t ; the maximum is indicated with a small arrow and labelled in bold. (Bottom row) Contours of magnetic flux 
t/j/tpi,, at the labelled times (solid), and at t = (dashed). The dotted contours meet their solid equivalents on the left vertical axis. Note 
that the scale of the r axis differs between Figures |(a)[)(d)| The change in the contour scale between Figure |(e)| and Figures |(f)jj(h)| is 
due to the difference in if) between t = 0, where the magnetic field is dipolar, and subsequent times when the magnetic field is distorted. 



distribution seen in Figure 7(a) This is not surprising. Mat- 
ter is fed in at r = r m i n with zero velocity, as discussed in 



section [2T1 It expands outward due to the pressure gradient 
created as matter piles up at the injection radius; since we 
are injecting quasistatically, the pressure gradient is small. 
On the other hand, the weight of the massive overburden 
(A/base ~ lOAia) presses down on the added material. The 
magnetic flux is displaced [Figure [8(e)| , but much less than 
for the hard mountain. Field lines remain tied to the inner 
boundary, bending away in its immediate vicinity (because 
the slug of injected matter does not rise). Above this layer, 
the field lines of the initial and final states remain largely 
parallel. 

The structure of the sunk mountain grown at r — R-,, 
[Figure |8(c) differs from the other two cases. The contour 
10, tracked above, starts at the pole, remains 
2ho below the injection radius, bends 



pX a /p H „ = 
virtually flat at 
sharply inward near the equator, moves directly toward the 
inner boundary, then curves back towards the pole, cross- 
ing it again at ~ 6ho below the injection radius. For the 
previous two mountains, grown from r — r m i u , the angular 
variation in density increases with altitude. Here the reverse 
is true: the angular variation density decreases with increas- 
ing r, up until r < R t , with the greatest variation within 
~ 4fto of the inner boundary. While the two sinking scenar- 
ios differ in their final distributions of accreted (as opposed 
to total) density, their final distributions of magnetic flux 
[Figures |8(e) | and |8(f )| are very similar. 

The mountain sinks three times further into the fluid 



base at the pole than at the equator. This is consistent with 
how mass is injected according to equation Q; the input 
flux is ~ 20 times greater at the pole than at the equator. In 
addition, the magnetic field guides accreted matter sideways 
as field lines flatten across the surface towards the equator, 
whereas matter at the pole can sink inward readily along 
almost vertical flux tubes [e.g. the contour t/j/i'* — 6 x 10~ 2 
in Figure [8(f)] . 

The density contours bunch together along the under- 
side of the mountain, spanning five orders of magnitude; 
the injected matter does not sink below this boundary. The 
lowest of the bunched density contours, pX a /p_R„ = 10~ 3 , 
never reaches the inner boundary; the mountain is floating 
in isostatic equilibrium with the surrounding fluid base. The 
contour rises to only ~ lho above R* at the equator and 
~ 3ho at the pole; in contrast it reaches R* — r< 6ho at the 



pole. The path of this contour, if overlaid on Figure 8(a) 
would trace densities between 2 x 10 -3 and 50pR^ . Finally, 
note that 0.3% of the mountain mass is above the stellar 
surface. Compared to the other two scenarios, the structure 
of <S(i?*, 10) is perhaps more reminiscent of an "iceberg". 

Ultimately we are interested in the final distribution of 
the total mass, that is, the accreted matter, pX a , plus the 
fluid base it displaces, p(l — X a ). Does injection at r m m or 
R* make a difference? Figure [9] displays the absolute, nor- 
malised difference [pr min — P-R*|/|Pr min + Pr„ \ between the 
total densities in the two sinking scenarios as a grayscale 
plot. The largest differences occur at r > i?*, where there is 
little mass, and are therefore unimportant. For r < R+, the 
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Figure 8. Comparison of the hydromagnetic structure of hard- and soft-surface mountains. Contours of accreted density pX a /pij t (top 
row) and magnetic flux ip/ipi (bottom row) of mountains 7i(10), iS(r m ; n , 10), and S(R+, 10), at time t = T a . Details are as for Figure[7] 
Note that the scale of the r axis differs between Figures |(a)))(c)| and between Figures |(d)[)(f )| 



difference peaks near the pole but remains less than ~ 0.8%. 
In other words, despite the difference in the final distribu- 
tion of pX a between the two injection scenarios (emphasised 
by the overlaid contours) , the final distribution of p is essen- 
tially the same. Injecting at r m i n or makes no difference, 
because the soft base readjusts in each case to yield the 
same overall equilibrium state. This is an important result. 
It confirms the robustness of the injection method and the 
argument presented in section [2~T1 in ideal MHD, the equi- 
librium state is independent of precisely where matter is 
initially injected. In practice, injection at r m i n seems prefer- 
able, because it reduces the simulation time (see Table [2} 
and improves numerical stability, as illustrated by the fail- 
ure of run <S(7?*,10 3 ). 



3.4 Magnetic line tying 

Finally, we investigate the assumption of magnetic line ty- 
ing. In Figure [TTJ1 we plot the angle (3(r,6) = sin -1 (Bg/\B\) 
between the magnetic field B and the radial unit vector as 
a function of 9 at r — r m in and r — R*. We also plot j3 
at the inne r boundary of a Grad-S hafranov mountain with 
M a = M c (|Pavne fc Melatosl [gPoj ) , and (3 for a dipole (in- 
dependent of radius). The S(R*) mountains give the same 
results as <S(r m i n ). 

The hard-surface mountains in Figure [lO (a) | behave like 
the Grad-Shafranov mountain at low M a but become in- 
creasingly dipolar as M a increases. This is expected; at low 
M a , the accreted mass stays close to the pole and distorts the 
magnetic field there. As Af a increases, the mountain spreads 



over a greater volume, and the magnetic field is distorted 
less at any particular point. The sign inversion close to the 
equator may be caused by numerical reconnection, or by the 
reflective boundary condition at 6 = n/2; further tests with 
a resisitive ideal-MHD solver are needed to make sure. 



The soft-surface mountains in Figure 10(b) are dipo- 
lar at r = Ri,, as expected, but at r = r m i n [Figure 10(c) 



they more closely resemble the Grad-Shafranov mountain. 
One might expect f3 to closely match a dipole at r = r m - m , 
given that the magnetic field lines are tied there, and we 
choose Mb a sc 3> M a in order to minimise sideways fluid 
displacements at the bottom of the soft base. It is unclear 
whether the magnetic distortions are artificial, because the 
injected slug matter enters from below and cannot expand 
upwards to match accretion from above (see section [5] for 
further discussion) . Alternatively, kinks in the magnetic field 
may be communicated rapidly down to arbitrary depths by 
Alfven waves, even though the Alfven speed oc p _1//2 de- 
creases rapidly with depth. If so, the high breaking st rain 
of the solid, conducting crust l|Horowitz fc K adau 200j|) as- 
sumes even greater importance in enforcing line tying. 

We argued, in section 12.11 that the final equilibrium 
state of the mountain is independent of ri n j. In general, a 
given total M a and injected mass flux dM & /dij) does not de- 
fine a unique ideal-MHD equilibrium. Matter injected from 
above spreads sideways faster than it sinks, like a layered 
cocktail drink, while a slug of matter injected from below 
forces the base sideways without much movement at the sur- 
face [compare Figures [8(c)| and |8(b)| . Conceivably, therefore, 
ZEUS-MP may converge on different equilibria depending 
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Figure 10. Angle between the magnetic field B and the radial 
unit vector, plotted versus colatitude 8 at |(a)| r = r m j n = i?*, |(b)| 
r = R*, and [(c)] r = r min , for Af a /M c = 1 (dotted), 10 
(dot-dashed), 10 2 (dashed), and 10 3 (solid). Results from this 
paper are plotted in b lack. Plotted in gray are /3(r, 8) for a 
iPavne fc M elatos (2004) Grad-Shafranov mountain (dotted), and 
for a dipole (solid). 
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Figure 11. Ellipticity as a function of time, for M a /M c = 1, 10, 
10 2 , 10 3 (top left to bottom right), and for mountains "H (solid), 
5(r m ; n ) (dashed), and <S(H*) (dotted). 



2.0 



on rinj. The results of section 13731 engender confidence that 
the mountain structure does not depend on ri n j; the issue 
is not definitively settled, however, for the following subtle 
reason Q 

Consider a polar field line in the bottom row of Fig- 
ure [Jj As accretion proceeds, it bends towards the equa- 
tor until it touches the corner (r, 8) = (r max,7r/2). At that 
point, it instantaneously snaps through some nonzero an- 
gle, from B r 7^ (free boundary at r = r max ) to B r — 
(reflecting boundary at 9 = 7i"/2). Effectively, this corre- 
sponds to a dissipative, reconnection-like event occurring 
just outside the simulation volume, artificially pinching off 
magnetic loops[f| Such a process is irreversible. Furthermore, 
it acts differently on the sequence of quasistatic equilibria 
that ZEUS-MP hypothetically passes through during slow 
accretion from above and below, because sideways spread- 
ing happens at different altitudes in the two cases. 

In the runs presented in this paper, the density in the 
vicinity of the corner (r max ,n/2) is tiny, as is the mass ef- 
flux through the boundary r = r max (see section I3.1|l . One 
can therefore argue plausibly that the irreversible dissipation 
at (r max , 7r/2), while it exists in principle, does not signifi- 
cantly affect the final state. There is a chance, however, that 
if one adds material slowly from above, reconnection (where 
numerical or real) pinches off one small magnetic loop af- 
ter another at the equator, as in the Earth's magnetotail . 
Resisitive MHD simulations by IVigelius fc Melatos! (|2009bl ) 
do not show such behaviour, but they mostly started from 



4 Sterl Phinney, private communication. 

5 The effect can be magnified in ZEUS-MP by increasing the 
cell size close to v — T'ma.x! eventually ZEUS-MP aborts when B r 
diverges close to the (r max ,7r/2) corner. 
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Figure 12. Ellipticity as a function of M a , at time t = T a , for 
mountains H (circles), iS(r m ; n ) (triangles), and <S(iJ*) (squares). 



preformed Grad-Shafranov equilibria instead of growing the 
mountain from scratch. A more careful consideration of this 
issue is required for future simulations. 



4 MASS QUADRUPOLE MOMENT 

The distorted hydromagnetic equilibria in section [3] have 
an associated mass quadrupole moment, with principal axis 
along the pre-accretion magnetic axis, which is quantified in 
terms of the ellipticity 



dr r 



r/2 



d9 sin 6>(3 cos 2 6>- l)p(t,r,8), (5) 



with I Z z — 2Af*i? 2 /5. Figure [TT] shows e with respect to time 
as the mountain grows. All mountains achieve a nonzero el- 
lipticity at t — T a , which decreases negligibly thereafter, 
confirming the mountains are stable. The time taken for 
the hard-surface mountains to converge to their equilibrium 
values of e decreases with M a , from ~ T a (M a = 1M C ) 
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Table 3. List of important physical parameters of accreting neu- 
tron stars (top part), and a summary of the results of the simu- 
lations presented in this paper (bottom part). 



Quantity 


Value/ Range 


Reference 


accreted mass 


10" 4 - O.8M 


1, 2, 4 


accretion timescale 


10 4 - 10 6 yr 


1, 7 


density of crust 


1Q 9 _ 1Q 14 g cm -3 


5, 10 


depth of crust 


~ 1000 m 


5, 10 


initial magnetic field 


10 12 - 10 13 G 


3, 8, 9 


temperature 


10 8 - 10 9 K 


6, 10 


cllipticity 


5 x 10~ 5 - 2 x 10" 4 


11 


effect of sinking 


e reduced by 25 - 60 % 


11 
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to ~ 0.1T a (M a = 10 3 Af c ); this is consistent with the de- 
creased confinement of the mountain by the magnetic field, 
i.e. its increased ability to spread. The ellipticities of the 
soft-surface mountains for the two injection scenarios are vir- 
tually identical; even the incomplete run <S(i?*,10 3 ) closely 
follows (S(r m i n , 10 3 ) up until failure. As in section 13.31 the 
final density distribution is independent of the injection pro- 
cedure. 

The origin of the uneven behaviour of the ellipticity of 
W(10 3 ) [Figure! ll(d)| * s unknown; we note, however, that its 
functional dependence on t is similar to that of M a [equa- 
tion and Figure [2J, and therefore it is likely that the 
rise and fall of the ellipticity is due to the reconfiguration 
of the 7i(10 3 ) mountain in response to the changing accre- 
tion rate. We note similar undulations in the kinetic energies 
[Figures |5(b)| and |5(c)| and, to a lesser degree, in the ellip- 



ticity of H(10) [Figure 11(b) 



Figure [12] shows e (in black) at t = T a as a function of 
M a . The ellipticity of the hard-surface mountains rises by 
a factor of ~ 2.5 as M a increases from 1M C to 10M C and 
flattens thereafter, rising by a further 5% as M a increases 
from 10M C to 10 3 M C . Importantly, accreting further matter 
does not reduce e; the mountain does not smooth itself out. 
The ellipticities of the soft-surface mountains rises by ~ 60% 
per decade in M a . 

Figure [12] clearly quantifies the effect of sinking: e de- 
creases, relative to the hard-surface scenario, by ~ 50% 
at M a = 1M C , ~ 60% at M a = 10Af o , and ~ 25% at 
M a = 10 3 M C . 



5 DISCUSSION 

In this paper, we simulate the growth of magnetically con- 
fined mountains on an accreting neutron star, with realistic 
masses M a < 0.12Mq, under the two scenarios where the 
mountain sits on a hard surface and sinks into a soft, fluid 
base. In the latter scenario, we confirm that the final equi- 
libriun state is independent of the altitude where matter is 
injected. We find that the ellipticity of a hard-surface moun- 



tain does not increase appreciably for M a > 10M C , saturat- 
ing at ~ 2 x 10~ 4 , whereas the ellipticity of a soft-surface 
mountain continues to increase from M a = 10A/ C to 10 3 M C . 
Sinking reduces the ellipticity by up to 60% relative to the 
har d-surface value. 

IChoudhuri fc Konarl l|2002l ) developed a kinematic 
model of accretion, which treats sinking in a different (but 
complementary) way to this paper. An axisymmetric mag- 
netic field is evolved under the influence of a prescribed ve- 
locity field, which models the flow of accreted matter from 
pole to equator, where it submerges and moves towards the 
core (see their Figure 1). Ohmic diffusion is included, but, 
for a subset of the results (where the resistivity r\ — 0.01), it 
is negligible, permi tting a direct compar is on wit h this paper. 

Figure 5 of IChoudhuri fc Konarl 1)20021 ) shows the 
evolved configuration of an initially dipolar field which per- 
meates the entire star. We compare to Figures 8(e) and 8(e) 



of this paper. In both models, the magnetic field is distorted 
significantly by accr eted matter spreading towards the equa- 
tor (r m < r < r s in Choudhuri & Konai; the entire simula- 
tion in this paper). In IChoudhuri fc Konarr s work, the mag- 
netic field is completely submerged beneath the surface and 
confined to the zone where the submerged accreted mat- 
ter flows back towards the pole. In this paper, magnetic 
field lines still penetrate th e surface, implying le ss effective 
screening. Within the core, I Choudhuri fc Konarl ' s magnetic 
field remains relatively undisturbed. Magnetic line-tying is 
not enforced, but the prescribed radial flow within the core 
naturally restricts the sideways displacement of the mag- 
netic field there. If there were sideways motion of the mat- 
ter within the core, it would modify the degree of mag- 
ne tic screening, but ne ither our simulations nor the results 
of IChoudhuri fc Konarl show evidence for such motion. Ex- 
tending our simulations deeper into the star to include the 
core and explore this possibility properly would be a tech- 
nical challenge; for instance, we would need to incorporate 
a more realistic equation of state and track even more dis- 
parate equilibrium time-scales. 

To explain the narrow rang e in the rotation frequen - 
cies of low-mass x-ray binaries l|Chakrabartv et al.l I2003T I. 
it is proposed that the stars radiate angular moment in 
gravita tional waves at a rate which b alances the accretion 
torque (Wagoner 1984; Bildstcn 1998). Magnetic mountains 
are one of a number of physical mechanis ms proposed for the 
associa ted permanent quadrupole; see IVigelius fc Melatosl 
l|2009al ) and references therein. The relationship between e 
and the rotation frequency / predicted by torque balance is 



-2/5 



Thus, the 25% to 60% reduction in e due to sink- 



ing calculated in this paper increases / by 12% to 44%, all 
other things being equal. This goes some way towards bring- 
ing magnetic mountain ellipticities down to a level consistent 
with the data, but there is still a long way to go. Observa- 
tions to date have found 45 Hz < / < 620 Hz for burst 
oscillation sources and 182 Hz < / < 598 Hz for accreting 
millisecond pulsars, implying 6.6 x 10~ 9 < e < 4.6 x 10~ 6 
and 7.2 < 10~ 9 < e < 1.4 x 10~ 7 respectively. Conversely, 
the ellipticities of sunk mountains calculated in this paper, 
3.5 x 10 -5 < e < 1.5 x 10~ 4 , imply 11 Hz < / < 20 Hz. 
Clearly, other relaxation mechanisms, like Ohmic diffusion, 
must also be playing an important role in reducing e, as the 
observed / require. 

The reduction in e by sinking also reduces the gravi- 
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tational wave strain (e.g. I Abbott et alj [20071 ) . h tx e/ 2 , by 
6% to 17%. This is unlikely, by itself, to rule out the detec- 
tion of gravitational waves from low-mass x-ray binaries by 
ground-based interferometric detectors; assuming the signal 
can be coherently integrated, the loss in h can be compen- 
sated for by an increase in the observation time oc h~ 2 of 
13% to 45%. Other difficulties associated with the detection 
of gravitational waves from low-mass x-ray binaries, such 
as poorly known orbital parameters and accretion-induced 
phase wandering (TWatts et al.ll2008h . are likely to be more 
important. 
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APPENDIX A: MATCHING A FLUID BASE TO 
A GRAD-SHAFRANOV MOUNTAIN IN 
ZEUS-MP 

We atte mpted to incorpor a te a fl uid b ase into the frame- 
work of iPavne fc Melatos! (|2007f ) and IVigelius fc Melatos! 
( 2008) in the following ad-hoc manner. Starting with a Grad- 
Shafranov equilibrium loaded into ZEUS-MP, we extend the 
inner simulation boundary, initially at r = R*, inwards to 
create a region r m i n < r < R+, containing the fluid base. 
The magnetic field B in this region is initialised to a dipole. 
At r = i?*, B r matches perfectly, but Be i s disc ontinu- 
ous [see, e.g., Figure 2 of IPavne fc Melatos! |2004)]. The 
initial density p(t — 0, r, 9) is chosen to match the Grad- 
Shafranov density pgs(t, 9) at r = i?*, and to match an 
isothermal, non-self-gravitating profile within r < R±. A 
number of ad-hoc choices of p(t = 0, r, 6) were tried, e.g. the 
maximum of pgs(R*,9) and p' Rt exp[GM*(r _1 — R^ )/cj], 
with p' Ri = mine pgs(R*, 9). When the combined Grad- 
Shafranov mountain and fluid base are evolved in ZEUS-MP, 
the results are undesirable. Except when i?* — r m i n <C ho, the 
fluid base is sufficiently far from equilibrium to completely 
disrupt the Grad-Shafranov mountain, which collapses over 
a short time-scale ~ to- 



APPENDIX B: LOGARITHMIC RADIAL GRID 
SPACING 

The logarithmic grid spacing in r is determined as follows. 
The N r + 1 radial cell boundaries r m i n = ro, n , r% , . . . , rjv r = 
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1 max 

are given by r n+ i — r„ + Ar n , where 

N r -1 

^ Ar„ = r max - r min , (Bl) 
Ar n+1 fAr Nr -i\ 



Ar„ \ Ar 

and ArN,,-i/Ara is the ratio of the maximum to minimum 
radial grid spacing. The values of Arjv r _i/Aro used in the 
simulations presented in this paper are given in Table [2] 



APPENDIX C: CUSTOM INJECTION 

We add a new subroutine to ZEUS-MP which is called at 
the beginning of each time-step St. Within the subroutine, 
the density p(t, r, 9), mountain concentration X a (t, r, 9), and 
velocity v(£, r, 9) of a grid cell within the injection region (at 
point (r, 9) with size Sr x 89) are updated, as follows: 

p(t + 6t, r, 9) = p(t, r, 9) + 5p(t, r, 9) , (CI) 

+ = p(t,^)x(t,^) + Mt, r , g) 

p(t + dt,r,&) 

+ v(t,r,9)[l-X4t + 5t,r,9)]. 
The density increment is given by 
M 

Sp(t,r,9) = ^-^l(t,r,9); (C4) 

the factor of 2n comes from the size of the grid cell in the <f> 
dimension. The function 

1 pH fr-2 re + SB 8 s M 

I{t ' r ' 9) = J7j tl dt J ri drr2 l Msine 8W^9^ 

(C5) 

integrates the injected flux given by equation (J2Jl; the con- 
stant 

(C6) 

ensures the correct normalisation. The times 

ti =min(i,T a ), (C7) 
t 2 = min(i + St, T a ) , (C8) 

give the intersection of the current time-step with the injec- 
tion time interval, and the radii 

n = min[max(r, r in j), r inj + 5r inj ] , (C9) 
r 2 = min[max(r + Sr, r in j), r inj + 5r inj ] , (CIO) 

give the intersection of the grid cell with the injection region. 



